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Abstract The Maxwell equations for the electromagnetic potential, supple- 
mented by the Lorenz gauge condition, are decoupled and solved exactly in de 
Sitter space-time studied in static spherical coordinates. There is no source 
besides the background. One component of the vector field is expressed, in its 
radial part, through the solution of a fourth-order ordinary differential equa- 
tion obeying given initial conditions. The other components of the vector field 
are then found by acting with lower-order differential operators on the solu- 
tion of the fourth-order equation (while the transverse part is decoupled and 
solved exactly from the beginning). The whole four- vector potential is even- 
tually expressed through hypergeometric functions and spherical harmonics. 
Its radial part is plotted for given choices of initial conditions. We have thus 
completely succeeded in solving the homogeneous vector wave equation for 
Maxwell theory in the Lorenz gauge when a de Sitter spacetime is consid- 
ered, which is relevant both for inflationary cosmology and gravitational wave 
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theory. The decoupling technique and analytic formulae and plots are com- 
pletely original. This is an important step towards solving exactly the tensor 
wave equation in de Sitter space-time, which has important applications to 
the theory of gravitational waves about curved backgrounds. 



1 Introduction 

It is by now well known that the problem of solving vector and tensor wave 
equations in curved spacetime, motivated by physical problems such as those 
occurring in gravitational wave theory and relativistic astrophysics, is in gen- 
eral a challenge even for the modern computational resources. Within this 
framework, a striking problem is the coupled nature of the set of hyperbolic 
equations one arrives at. For example, on using the Maxwell action functional 



F ab F ab ^ d A x 



A I 



jointly with the Lorenz \T\ gauge condition 



V b A h = 0, 



(1) 



(2) 



one gets, in vacuum, the coupled equations for the electromagnetic potential 



(3) 



It was necessary to wait until the mid-seventies to obtain a major break- 
through in the solution of coupled hyperbolic equations such as (3), thanks 
to the work of Cohen and Kegeles [2] , who reduced the problem to the task 
of finding solutions of a complex scalar equation. Even on considering spe- 
cific backgrounds such as de Sitter spacetime, only the Green functions of 
the wave operator have been obtained explicitly so far [3] , to the best of our 
knowledge. 

Thus, in a recent paper [4], we have studied the vector and tensor wave 
equations in de Sitter space-time with static spherical coordinates, so that 
the line element reads as 



ds 2 = -f di 2 + r 1 dr 2 



'■(de 2 



sm 



(4) 



where / = 1—Hr, and H is the Hubble constant related to the cosmological 
constant A by H 2 — A, The vector field X solving the vector wave equation 
can be expanded in spherical harmonics according to [5] 



X = Y lm (0)e^-™» f (r)dt + h{r)dr 

I —i{u)t—m4 



'" r f2(r)Y lm (e) + f 3 (r) d) ' 1 



-i{ujt—m4 



sin(9" 



-r sin 0/ 2 (r) 



dO 



de 



dYi m 

de 



mf 3 (r)Y lm (0) 



(5) 
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where Yi m (0) is the ^-dependent part of the spherical harmonics Yi m (9, <j>). 
As we have shown in [4] , the function f 2 is decoupled and obeys a differential 
equation solved by a combination of hypergeometric functions, i.e. 



mo = r 1(2,2 

where 



U x r l F ^a_, a+; ~ + I; H 2 r 2 ^j +U 2 r- l - 1 F ^a+, a,-;~-l; H 2 r 2 

(6) 

U= -^r, a± = --(2ifi-3-2/±l). (7) 

At this stage, however, the problem remained of solving explicitly also 
for /o(r), /i(r), fz{f) in the expansion (5). For this purpose, Sec. II derives 
the decoupling procedure for such modes in de Sitter, and Sec. Ill writes 
explicitly the decoupled equations. Section IV solves explicitly for /o,/i,/3 
in terms of hypergeometric functions, while Sec. V plots such solutions for 
suitable initial conditions. Relevant details are presented in the Appendices. 
It now remains to be seen whether a technique similar to Sees. II and III can 
be used to solve completely also the tensor wave equation obtained in [4j. 
This technical step would have far reaching consequences for the theory of 
gravitational waves in cosmological backgrounds, as is stressed in [3], and we 
hope to be able to perform it in a separate paper. 



2 Coupled modes 

Unlike f 2 , the functions /o,/i and fa obey instead a coupled set, given by 
Eqs. (54), (55), (57) of [4j, which are here written, more conveniently, in 
matrix form as (our L = 1(1 + 1), and we set e = 1 in the Eqs. of 4 , which 
corresponds to studying the vector wave equation (3)) 



having defined 





A 3 \ //o\ 

P u r- 2 T X LC z /i = 0, 
C3 -P33 / \ fs J 


(8) 




2if2H 3 r 

A '= j • 


(9) 




r 


(10) 




d 2 d 

Poo = -7-5- + Qi~r + O2, 

dr z dr 


(11) 




d 2 d 
dr A dr 


(12) 




d 2 d 

P33 = -7-9 + Q5-7- + Q6, 
dr dr 


(13) 




Qi = C3 = -) 

r 


(14) 
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Q 2 H 2 L 



(15) 



q 3 = Ui-1±\, m 



r 



3/ 



UH 2 + v 7 -, (17) 



/2 y r 2 / / 

Q 5 = * (l - ~) , (18) 

Q 2 H 2 LI 

Q* = -p--^j, (19) 

where we have corrected misprints on the right-hand side of Eqs. (54) and 
(55) of [1]. With our notation, the three equations resulting from (8) can be 
written as 

(20) 

^A. ( 21 ) 

(22) 



Poo/o — 


-A 3 h, 


Pll/l = 


A 3 




P33/3 = 


-c 3 h. 



3 Decoupled equations 

We now express /1 from Eq. (20) and we insert it into Eq. (21), i.e. 

j nn / n I — 

A3 



— — Poo/o I - --p /o - 72^/3- 



(23) 



Next, we exploit the Lorenz gauge condition (2), i.e. [I] 



1 



Lfa = r z f— I -^-Poo/o ) - 2r(l - 2/) ( -—Poo/o ) + i 



1 



A, 



/ 



/o, (24) 



and from Eqs. (23) and (24) we obtain, on defining the new independent 
variable x — rH, the following fourth-order equation for /q: 



£ + + P 2 (x)^ + 2fc0«0| + Po(x) 



where 



/o(x) = 0, (25) 



(26) 



x A {x 2 - l) 4 

b (x) = L{L-2) + 2L{2-L-Q 2 )x 2 +\Q i + 'iQ 2 + L(L + 2([2 2 ~l))\x i , (27) 



D . , A{ fi 2 + L - 2 + Qx 2 ) 
B1{X) ^ xjx^lj 2 : 



(28) 
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B 2 (x) 



-L + (f2 2 + L- U)x 2 



18x 4 



x 2 (x 2 — l) 2 
x{x z — 1) 



(29) 



(30) 



Eventually, fi and F 3 = Hf 3 can be obtained from Eqs. (20) and (24), i.e. 



AW 



i (1 - x 2 ) ( d 



2Q 



\'dx 2 



fi 2 



L 



2 d 



x dx (1 — x 2 ) 2 x 2 (\ — x 2 ) 



LF 3 (x) 



x 2 (l~x 2 )^-2x(2x 2 -1) 
dx 



fi(x)+W- 



fo(x), (31) 



:fo(x). (32) 



(1 - a; 2 )' 

Our fi and f 3 are purely imaginary, which means we are eventually going to 
take their imaginary part only. Moreover, as a consistency check, Eqs. (31) 
and (32) have been found to agree with Eq. (22), i.e. (22) is then identically 
satisfied. 



4 Exact solutions 

Equation (25) has four linearly independent integrals, so that its general solu- 
tion involves four coefficients of linear combination C\, C%, C 3 , C4, according 
to (hereafter, F is the hypergeometric function already used in (6)) 

f Q (x) = c x x+\i - x 2 y^F (- l ~n - l -, - l -n + \ - i; \ - z ; r 

+ c 2X - x -\i - x 2 )-^ 2 f (- l -n + 1 - -, —n - - - -■ - - l; r 
v ; V 2 2 2 2 2 2 

+ c 3 x\i x 2 )-^f ^-n + i-'-n + l + l-l + i.x 2 

+ C,x\l - x 2 y¥i F {-\0 + l + l- l -n + \ + l --\ + l-, x 2 ^j .(33) 

Regularity at the origin [x = should be included, and we recall that the 
event horizon for an observer situated at x = is given by x = 1 [6]) implies 
that Ci = Ci =0, and hence, on defining 

i ^ I , *^3Z, 3, 

ai = --n + -, bi = --n+- + -, d ± = - + i, (34) 

we now re-express the regular solution in the form (the points x = 0, 1 being 
regular singular points of the equation (25) satisfied by /o) 

fo{x) = x l (l-x 2 )-i n [c 3 F( ai , h; df, x 2 )+dF( ai + l, h-l; d l5 x 2 )] , (35) 

where the second term on the right-hand side of (35) can be obtained from 
the first through the replacements 

C 3 — > C4,, ai — > ai + 1, 61 — > 61 - 1, 
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and the series expressing the two hypergeometric functions are conditionally 
convergent, because they satisfy Re(c — a — b) — if2, with 

a = Oi,oi + l; 6 = &i, &i — 1; c=d\. 

Last, we exploit the identity 



-^-Fta, b; c; z) = —F(a + 1, b + 1: c+ 1; z) 
dz c 



(36) 



to find, in the formula (31) for f\(x) 
d 



dx 
2axb 



f (x) = C 3 ix c - 1 (l-x") 



17-1 



1{1 - x 2 ) + if2x 2 F(a u bx;di;x 2 ) 



V ^x l+1 {\ - x 2 )-i n F{ ai + 1, h + 1; d 1 + 1; x 2 ) J 
+ {C 3 -» C 4 , ax -> oi + 1, h -» 6i - l}. 



(37) 



It is then straightforward, although tedious, to obtain the second derivative 
of /o (see Eq. (38) of Appendix A) in the equation for fi, and the third 
derivative of fo in the formula (32) for Hfe. The results are exploited to plot 
the solutions in Sec. V. 

In general, for given initial conditions at a G [0, 1[, one can evaluate C3 
and C4 from 

f (a) = 13, f Q (a) = 7, 
i.e. C r 3 = C 3 (/3,7),Ci = C 4 C8,7). 



5 Plot of the solutions 

To plot the solutions, we begin with /o as given by (35), which is real-valued 
despite the many i factors occurring therein. Figures 1 to 3 describe the 
solutions for the two choices C3 = 0, C4 = 1 or the other way around, and 
various values of I and fl. 

We next plot fi/i and F$/i = Hfe/i by relying upon (31) and (32). As 
far as we can see, all solutions blow up at the event horizon, corresponding 
to x = 1, since there are no static solutions of the wave equation which are 
regular inside and on the event horizon other than the constant one [6] . 



A Derivatives of fo 

The higher-order derivatives of fo in sections III and IV get increasingly cumber- 
some, but for completeness we write hereafter the result for fo(x), i.e. 

d^ Mx) = C - A { MdTTT) X {1 ~ X) n«i+2,6i + 2;dx + 2;x) 

di 
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Fig. 2 Regular solution (35) for f with C 3 = 1, C 4 = 0, 1 = 2, . . . 6 with 1? = 2 
(left figure) and J? = 4 (right figure). Increasing values of I correspond to more 
peaked curves on the right part of the plots. 
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0.2 0.4 0.6 0.8 1 



Fig. 3 Regular solution (35) for f with C 3 = 0, C 4 = 1, 1 = 10 and Q = 10. 




Fig. 4 Regular solution (31) for fi/i with C3 = \,d = (left figure) and C3 
0, C 4 = 1 (right figure) for 1 = 2 ... 6 and Q = 2. 



x [(2Z + 1)(1 - a; 2 ) + 2iflx 2 ]F{a 1 + 1, 61 + 1; di + 1; a; 2 ) 



+ x ; - 2 (l-^ 2 )-^- 2 



Z(Z - l)(x 2 - l) 2 - ^(a: 2 - l)(lx + 2(1 + l)x 2 ) 



+ (2iQ - f2 2 )x" 



F(a u bi 



;di;a; 2 )| 
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Fig. 5 Regular solution (32) for F$/i with C3 = 0, C4 = 1 (left figure) and C3 
1, Ca = (right figure) for / = 2 ... 6 and « = 2. 



+ < C3 — > C4, 01 — ► ai + 1, 61 — » 61 — 1 



(38) 



B Special cases: I — and £ = 1 

We list here the main equations for I — 0, f for completeness. These equations do 
not add too much to the above discussion and hence we have decided to include 
them in the appendix. 



B.l The case I = 0,m = 

In this case we have Y(6) = (47r) -1 ^ 2 = constant and the only surviving functions 
are /o and /1 . The main equations in [4] reduce then to (recalling that our Eq. (3) 
requires setting e = 1 in [4]) 



d 2 f _ 2 dfo 



Q 2 , 2iQx , 
JO + TZ7> TT/i 



dx 2 x dx (x 2 — l) 2 JU ' (x 2 — 1) 
and the Lorenz gauge condition (24) which now becomes 



(39) 



f2fo = i(x — 1)' 



dfx 2i(x 2 - l)(2x 2 - 1) 



dx 



+ 



h 



(40) 



These equations can be easily separated and explicitly solved in terms of hyperge- 
ometric functions. The details are not very illuminating and are therefore omitted. 
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B.2 The case I = l,m = 0, 1 
In this case we have 

Y = \ — cose l = \,m = Q 

V 47T 

Y = - l = l,m = l. (41) 

However, by virtue of the spherical symmetry of the background, the equations for 
both cases I = 1, m = and I = 1, m = 1 do coincide. We have 

d 2 / _ _2^/o _ (f? 2 x 2 + 2x 2 -2) + 2iQx 
da; 2 x dx a; 2 (a; 2 — l) 2 (x 2 — 1) 

d 2 /i _ 2(3x 2 - 1) dfi (0 2 x 2 + 4x 4 - 4) UQx 



dx 2 x(x 2 - 1) dx x 2 {x 2 - l) 2 ^ + (x 2 - 1)3 *° 

4F 3 
x 3 (x 2 — 1) ' 

d 2 / 2 _ 2(2x 2 - 1) d/ 2 _ (J? 2 x 2 + 2x 4 -2) 
dx 2 x(x 2 — 1) dx x 2 (x 2 — l) 2 

d 2 F 3 _ 2x dF 3 (fi 2 x 2 + 2x 2 - 2) 2 



dx 2 (x 2 -l) dx x 2 (x 2 -l) 2 Fs x fl ' (42) 

To this set one has to add the Lorenz gauge condition (24), which now reads 

Ofo = l[x2 - 1)2 1: + t^^ 3 + /ia;(2a;2 " 1)] • (43) 

Once more, the detailed discussion of this case can be performed by repeating 
exactly the same steps as in the general case, and is hence omitted. 
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